Is there a 4.5 PeV neutron line in the cosmic ray spectrum? 
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Recently we presented a model to fit the cosmic ray spec- 
trum using the hypothesis that the electron neutrino is a 
tachyon. The model predicted the existence of a neutron flux 
in the cosmic rays in a narrow region centered on E = 4.5±2.2 
PeV. The published literature on Cygnus X-3 reveals just such 
a 6a spike of neutral particles centered on E = 4.5 PeV. A 
second prediction of the model concerning integrated neu- 
tron fluxes at several energies also is consistent with published 
data. A specific further test of the model is proposed. 

PACS: 14.60.St, 14.60.Pq, 95.85.Ry, 96.40.De 



I. INTRODUCTION 

Recently we presented a model to fit the high energy 
cosmic ray spectrum using the hypothesis that the elec- 
tron neutrino is a tachyon. p] A good fit to the spectrum 
was obtained using |m„| = \J — m 2 = 0.5 ± 0.25 eV/c 2 . 
The signature prediction of the model is the existence 
of a neutron flux 'spike' in the cosmic rays centered on 
E = 4.5 ± 2.2 PeV, and having a width AlogS = 0.1 
(FWHM) . Although the existence of neutral cosmic rays 
from point sources remains a highly controversial subject, 
we report here that an examination of the published lit- 
erature on cosmic rays from Cygnus X-3 reveals just such 
a hitherto unreported neutral particle spike centered on 
E = 4.5 PeV with a level of statistical significance of 6a. 
An additional prediction of the model that the integrated 
flux of neutrons above 0.5 EeV should be 0.048 percent 
that above 2 PeV is also consistent with results from two 
out of three experiments. 

Although few physicists have taken tachyons seriously 
since they were first proposed in 1962 (2|, their exis- 
tence is clearly an experimental question. In 1985 Cho- 
dos, Hauser and Kostelecky ||, suggested that neutrinos 
were tachyons - an idea that is consistent with experi- 
ments used to determine the neutrino mass. Chodos et 
al. J|,|5| also suggested a remarkable empirical test of the 
tachyonic neutrino hypothesis, namely that stable parti- 
cles should decay when they travel with sufficiently high 
energies. Consider, for example, the energetically for- 
bidden decay p — > n + e + + v e . In order to conserve 
energy in the CM frame the neutrino would need to 
have E < 0. But tachyons with m 2 < have E < p, 
and therefore the sign of their energy in the lab frame 
Eiab = l{E + (3pcos9) will be positive for a proton ve- 
locity (3 > (3 t h = —E/pcos0. With the aid of a little 



kinematics it can easily be shown that the threshold en- 
ergy for proton decay is E t h ~ 1.7|m^| _1 PeV, with \m v \ 
in eV. 

Thus, if neutrinos are tachyons, energetically forbidden 
decays become allowed when the parent particle has suffi- 
cient energy - in seeming contradiction with the principle 
of relativity that whether or not a process occurs should 
not depend on the observer's reference frame. That con- 
tradiction is only an apparent one, however, because 
what appears to the lab observer as a proton decay emit- 
ting a neutrino appears to the CM observer as a proton 
absorbing an antineutrino from a background sea. 



II. COSMIC RAYS 

Since cosmic rays bombard the Earth with energies far 
in excess of what can be achieved in present day accel- 
erators, it is natural to ask whether any evidence for a 
process such as proton decay exists there at very high 
energies. One striking feature of the cosmic ray spec- 
trum is the "knee" or change in power law that occurs at 
E « 4 PeV. Various two-source mechanisms have been 
suggested to account for this spectral feature, but some 
researchers have identified it as arising from a single type 
of source. H In 1992 Kostelecky [Q] suggested that for a 
tachyonic neutrino mass |m„| Q.ieV, the proton decay 
threshold energy occurs at the knee of the cosmic ray 
spectrum, and could explain its existence. The idea is 
that cosmic ray nucleons on their way to Earth would lose 
energy through a chain of decays p—>n^p—>n— >■■■■, 
which would deplete the spectrum at energies above E t h- 
However, Kostelecky regarded the existence of the knee 
by itself as insufficient evidence for the tachyonic neutrino 
hypothesis in view of other more conventional explana- 
tions of the knee of the cosmic ray spectrum. He also 
did not attempt to model the spectrum, nor mention the 
signature neutron spike. 

Recently this author has developed a tachyonic neu- 
trino model that fits a number of features of the cosmic 
ray spectrum in addition to the knee. |l| These include 
the existence and position of the "ankle" (another change 
in power law at E ss 6 EeV), the specific changes in power 
law at the knee and ankle, the changes in composition of 
cosmic rays with energy, and the ability of cosmic rays 
to reach us above the conjectured GZK "cutoff." 
Although the fit to the cosmic ray spectrum was a good 
one, the model is highly speculative, because it is at vari- 
ance with conventional wisdom about cosmic rays and it 
arbitrarily assumed that the decay rate for protons (for 
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E > E t h) was far greater than that for neutrons. 

Nevertheless, the model did make the striking predic- 
tion of a cosmic ray neutron flux in a narrow range of 
energies just above E t h — a neutron "spike." The pile up 
of neutrons in a narrow interval just above E t h is a con- 
sequence of the fractional energy loss of the nucleon in 
proton decay becoming progressively smaller, the closer 
the proton energy gets to E t h- The position of the pre- 
dicted cosmic ray neutron spike depends on the value as- 
sumed for \m v \. From the fit to the cosmic ray spectrum 
we found |m„| = 0.5 ± 0.25 eV/c 2 , and hence we pre- 
dicted a neutron spike at E = 4.5 ± 2.2 PeV. In fact the 
model predicted that most nucleons should be neutrons 
for E > E t h, because it was assumed that as nucleons 
lose energy in the p — ► n — ► p — ► • • • decay chain, the 
lifetime and hence the decay mean free path for neutrons 
is far greater than for protons, and so nucleons above 
Eth would spend nearly all of their time en route as neu- 
trons. But, the model also predicts that for energies 
above the spike the neutron component does not become 
an appreciable fraction of the total cosmic ray flux until 
around 1 EeV. While neutrons might reach Earth at EeV 
energies in conventional cosmic ray models, it would be 
difficult to understand any sizable neutron component at 
energies as low as E=4.5 PeV, where the neutron mean 
free path before decay would be only about 100 ly. In 
the present model, however, A = 1 cosmic rays can travel 
very many neutron decay lengths and still arrive as neu- 
trons because many steps of the p — > n — > p — > • ■ ■ decay 
chain occur for nucleons having energies above E t h- 



III. CYGNUS X-3 DATA 

One way to look for a neutron flux would be to find 
a cosmic ray signal that points back to a specific source, 
since neutrons are unaffected by galactic magnetic fields. 
Starting in 1983 a number of cosmic ray groups did, in 
fact, report seeing signals in the PeV range from Her- 
cules X-l and Cygnus X-3. At the time these signals 
were believed to be either gamma rays or some hitherto 
unknown long-lived neutral particle, since neutrons, as 
already noted, should not live long enough to reach Earth 
(except in the present model). Some of the experiments 
coupled detection of extensive air showers with detec- 
tion of underground muons. |Tl"| , |l2| ] The observed high 
muon intensity was found to be consistent with hadrons 
but not with showers induced by gamma rays. [^2|,^3| It 
was widely believed that the mass of the neutral parti- 
cle was m 1 GeV/c 2 . juj Thus, all the observed or 
conjectured properties of these particles were consistent 
with neutrons: neutral strongly interacting particles with 
to « 1 GeV/c 2 . 

Following a period of excitement in the 1980's, many 
researchers began to look critically at some of the ob- 
servations of ultra-high energy cosmic rays from point 
sources. This skepticism was based in part on the in- 



consistencies between results reported in different exper- 
iments. As Chardin and Gerbier have noted fl5|j , a num- 
ber of papers used data selection procedures that made 
direct comparisons difficult, e.g., using different phase 
intervals to make cuts, variously reporting the total flux 
or only the flux in a particular phase bin, and reporting 
only "muon-poor" events. Also, some papers appeared 
to inflate the statistical significance of their results. 

But, the most serious challenge to the idea of neutral 
particles in the PeV range from Cygnus X-3 and other 
point sources came from a trio of high sensitivity experi- 
ments [jl6|-|l8| that reported seeing no signals from point 
sources claimed earlier. In the most sensitive experiment 
of the three, the upper limit on the flux of neutral par- 
ticles from Cygnus X-3 above 1.175 PeV was far below 
the fluxes reported by those experiments claiming sig- 
nals earlier. [[19| There seems to be only two possibilities: 
either all the earlier experiments claiming signals were 
in error, or Cygnus X-3 and other reported sources all 
had turned off about the time improved instrumentation 
became available. Table I offers some support for the 
latter possibility, because (a) the phases of the signals 
are in rough agreement in three experiments, and (b) the 
integrated flux above a PeV does appear to systemati- 
cally decrease over time taking all experiments together. 
(Among those claiming signals only those claiming more 
than 4cr have been listed, and among those citing upper 
limits only those giving upper limits on the flux above a 
PeV have been listed.) The suggestion that signals from 
Cygnus X-3 have fallen with time was first raised by N. 
C. Rana et al. based on X-ray and gamma ray data in 
four different wavelength regions. [Q ln what follows, 
we make the "optimistic" assumption that earlier exper- 
iments were seeing real signals, and we consider to what 
extent those reports of signals from Cygnus X-3 support 
the prediction of a 4.5 PeV neutron spike. 

In the 1980's there were eight cosmic ray groups that 
cited fluxes in the PeV range of signals pointing back to 
Cygnus X-3, (some which were inconsistent as mentioned 
earlier.) In nearly all cases limited statistics required 
reporting the flux integrated over energy in only one or 
at most two energy intervals. 



Ref Years E in PeV 



[ 23| 
[221 
[M| 

[LSI 
[L9| 



Flux Stat. sig. 



Phase 



76-79 

78- 81 

79- 82 
86-88 

89 
90-95 



> 2 

> 1 

> 3 

> 1 

> 1 
> 1.175 



7.4 ±3.2 

< 3 

1.5 ±0.3 
2.7 ±0.5 

< 23 

< 0.1 



4.4a 



0.1-0.3 



5a 0.225-0.25 
4.7a 0.25-0.30 



TABLE I. Experiments reporting integrated fluxes (or up- 
per limits) in units of x 10~ 14 particles cm -2 sec -1 for Cygnus 
X-3 for PeV energies. Only experiments reporting nonspo- 
radic signals claimed to be at a level of more than 4a have 
been listed. The van der Klis and Bonnet-Bidaud ephemeris 
has been used for finding the phase interval in each case. 
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One group (Lloyd Evans et al. however, had good 
enough statistics to report fluxes in eight energy bins 
spanning the location of the predicted 4.5 PeV neutron 
spike, and it had an energy acceptance threshold near Eq 
— 1 PeV, which could give one energy bin before the spike 
itself. The signal seen by Lloyd Evans et al. from Cygnus 
X-3 did not appear until the data is selected on the basis 
of orbital phase determined from the X-ray binary's 4.79 
h orbital period, and the time of signal arrival. Lloyd- 
Evans et al. found that if they looked at the number 
of counts in 40 phase bins, one of these bins showed a 
sizable excess (73 counts when the average was 39). The 
information in Table II is taken from Lloyd-Evans et al. 
p2| , with the last column added by this author. Fig. 1 
displays the data in that last column. Wc would expect 
a flat distribution on the basis of chance, assuming that 
the signal were just a statistical fluctuation. In fact, av- 
eraged over all phases, the distribution must be flat and 
zero height, regardless of whether the signal is real or 
not. Note, that a spike appears centered on the value 
predicted by the tachyonic neutrino model, and that all 
the remaining bins have a flux consistent with zero. The 
gaussian curve drawn with arbitrary height in the fig- 
ure shows what would be predicted by the model given a 
neutron spike of width AlogE 1 = O.l(FWHM) and a 50 
percent energy resolution (Alog-E = ±0.176). Accord- 
ing to Lloyd-Evans, the actual resolution was probably 
around 50 percent, and very likely less than 100 percent 
p5[ We estimate the statistical significance of this spike 
occurring by chance by dividing the excess number of 
events in the two bins straddling 5 PeV by the square 
root of the expected number of events in those two bins: 
28.4/V22.6 = 6.0cr. It is interesting that in their arti- 
cle, Lloyd-Evans et al. displayed only the integrated flux 
/(> E) versus energy, and hence failed to mention the 
spike. Instead, they simply noted that the integrated 
spectrum appeared to steepen right after 10 PeV. 

How can we be sure that the spike seen in Lloyd-Evans 
et al. data is not an artifact of the data analysis or a sta- 
tistical fluctuation? Six standard deviations may seem 
interesting, but the original peak in their phase plot was 
far less impressive, particularly allowing for a "trials fac- 
tor" of 40, since such a peak might have been seen in 
any one of the 40 phase bins. Suppose that in fact the 
original peak in the phase plot were a statistical fluctua- 
tion, how could one then get a 6cr peak in the flux versus 
energy distribution for events in a specific phase bin? 
Clearly, such a peak would require some correlation be- 
tween energy and phase. This could in principle occur, 
because observed cosmic ray energy is correlated with 
declination angle, and hence with time of day. However, 
all cosmic rays in a given phase bin arrive at one of five, 
i.e., 24/4.79, times throughout the day, and those arrival 
times slowly advance from day to day, since the Cygnus 
X-3 period is not exactly divisible into 24 hours. Thus, 
over the years of data-taking each phase bin would sam- 
ple times of the day with an almost uniform distribution, 
making it difficult to see how a phase-energy correlation 



could occur. 



E (in PeV) Observed Expected Excess ±V Expected 

L3 16 TX9 2.1± 3.7 

3-5 34 16.4 17.6± 4.0 

5-11 17 6.2 10.8± 2.5 

11-18 4 2.4 1.6± 1.6 

18-36 3 4.3 -1.3± 2.1 

36-72 6 3.4 2.6± 1.8 

72-140 2 0.8 1.2± 0.9 

>140 0.6 -0.6± 0.8 



TABLE II. Observed and expected event counts reported 
by Lloyd-Evans et al. in differential energy bins for the phase 
interval 0.225-0.250. The last column has been added by the 
author. The "Expected" counts for each energy interval are 
based on the average over all phases. 
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0.5 1 1.5 2 

Log(Energy in PeV) 

FIG. 1. Data points from the last column of Table II plot- 
ted at the middle of each interval in log E in PeV. The gaus- 
sian curve centered on 4.5 PeV is what one would expect 
to find in Lloyd-Evans data, given a neutron spike of width 
Alog_E = O.l(FWHM), and a 50 percent energy resolution 
(A log.E = ±0.176) 



(It could be that at their source the phase and energy 
of cosmic rays are correlated, but in that case we would 
be dealing with a real source, not a statistical fluctuation, 
as hypothesized above.) 

Ideally, one would want to combine the Lloyd-Evans et 
al. data with that of other experiments in the PeV region 
to see if the spike either is destroyed or enhanced. Several 
problems arise with the other existing data, in which a 
signal is claimed from Cygnus X-3: one experiment used 
only "muon-poor" events |2q] , two experiments reported 
only the integral flux above some energy (no energy bin 
defined) [P7y2q| ; two reported the flux in an energy bin 
three times the width used by Lloyd-Evans 0,^91, and 
none was contemporaneous with Lloyd-Evans, thereby 
severely diminishing their utility. 

Aside from the spike, one other prediction of the tachy- 
onic neutrino model is that neutrons should also be seen 
as a significant and rising fraction of the cosmic ray flux 
above around 1.0 EeV. In fact, two cosmic ray groups 
have reported seeing neutral particles from Cygnus X-3 
having energies above 0.5 EeV with fluxes of 1.8 ± 0.7 
p(i| |, and 2.0 ± 0.6 ]3l|] , while a third group reporting 
merely an upper limit to the flux < 0.4 |32| - all in units 
of 10~ 17 particles cm~ 2 s . 

These measured fluxes above 0.5 EeV can be compared 
directly with the neutron flux predictions from the tachy- 
onic neutrino model, |IJ As noted previously, the ratio of 
the integral flux of neutrons above 0.5 EeV to that above 
2 PeV is predicted to be R = 4.8 x 10~ 4 . The predicted 
neutron flux for E > 0.5 EeV is then R times the mea- 
sured flux reported by Lloyd-Evans et al. for E > 2 PeV, 
or: R x 7.4 ± 3.2 x 10~ 14 = 3.5 ± 1.5 x 10~ 17 particles 
cm~ 2 s _1 , which is in quite good agreement with the two 
groups that measured a flux, rather than an upper limit. 
Although subsequent data accummulation by these two 
groups failed to show a signal from Cygnus X-3 |33l , that 



only adds additional support to the hypothesis that the 
source faded over time. 

If it is true that Cygnus X-3 and other point sources 
were active in the early 1980's and subsequently have 
turned off, is there any way to check whether there really 
is a 4.5 PeV neutron spike without waiting for specific 
sources to come back on? Without knowing where the 
sources are, the model can make no prediction of the 
anisotropy or the the angular distribution of sources of 
high energy cosmic rays. However, recall that the model 
predicts that all the cosmic rays include a 4.5 PeV neu- 
tron spike, not just those pointing back to the handful 
of possible sources looked at so far. Thus, if one selects 
events in a narrow energy band centered on 4.5 PeV, one 
could look at their arrival directions on the two dimen- 
sional map of the sky, and see if there is a noticeable 
clustering of points, which would indicate neutral par- 
ticles coming from specific sources. Moreover, if those 
sources were episodic, one should observe a nonuniform 
distribution in arrival times for events for a given source. 

Consider a specific example. The integrated flux in the 
4.5 PeV spike is 0.1 neutrons per m 2 -sr-s, which would 
give around 3 million counts over 5 years for an array of 
area 250,000 m 2 . If the array had an energy resolution of 
100 percent, it would also record a background count rate 
roughly four times as great in the energy bin centered on 
4.5 PeV. Suppose the angular resolution were A8 = 0.01 
rad, which would allow up to 4/A0 2 = 4 x 10 4 solid angle 
bins to be defined. Each bin would then have on the av- 
erage 400 background counts. Further suppose that the 
cosmic rays reaching Earth came from N point sources, 
then those solid angle bins pointing back to sources would 
have an average signal to background ratio: 10 4 /A^. Iden- 
tification of sources should then be possible, unless N 
were larger than the number of solid angle bins, and no 
subset of sources were appreciably brighter than others. 



IV. SUMMARY 

In summary, a highly speculative tachyonic neutrino 
model |IJ, which fits the cosmic ray spectrum well, pre- 
dicts a spike of neutrons at an energy where, given 
the neutron lifetime and distance to likely sources, very 
few should appear. A search through the literature for 
sources of neutral cosmic rays has identified a particular 
experiment with a favorable energy acceptance threshold, 
good enough statistics, and enough energy bins spanning 
the region of the neutron spike to test the prediction. The 
data do show a 6a spike located right at the predicted en- 
ergy, which was not identified in the original work. The 
failure of other subsequent more sensitive experiments to 
see a signal from Cygnus X-3 would seem to require that 
this source has since turned off - a possibility given some 
support by both time trends of data from different ex- 
periments, and data within the same experiments. The 
characteristics of the neutral particles from Cygnus X-3 
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seem to be consistent with neutrons rather than gamma 
rays, based on muon data from various experiments. For 
the EeV region, where the model also predicts neutrons 
(though not a spike) , two out of three experiments show 
a positive signal from Cygnus X-3, and they report a flux 
whose magnitude (relative to the flux in the spike) is well- 
predicted by the model. The hypothesis that the electron 
neutrino is a tachyon would seem to be supported, and 
it can be further tested without waiting for specific point 
sources to come back on. 
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